A  negative  rotation  of  the  vector  (or  positive  rotation  of  axes)  produces  the 
transpose  of  the  siatrix.  This  can  be  seen  directly  since  pairs  of  q's  in  the  off 
diagonal  terms  which  do  not  contain  q^  are  even  in  0  while  those  containing  q 
are  odd.  ° 

A  triad  of  orthogonal  vectors  can  be  introduced  by 


security  classification  of  this  PAGE  flWiwi  Pat*  gnur.dj 

REPORT  DOCUMENTATION  PAGE  before0 tompl  etwg^form 

1.  REPORT  NUMBER  JT  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 

TR  79-434  / 


«.  TITLE  (and  Subtitle) 

_gIEJK§.^TH0D  IN  LIGHT  OF  GENERALIZED  JULER 

■T&RAME*lfkS  0  "2S-  sS'  y 


of  report  ft  Ftmeo  covered 


Final  ^ 


17.  AUTMQRf+fr— 


6.  PERFORMING  ORG.  REPORT  NUMBER 


•.  contract  or  grant  numbert*) 


R.  N .j 


DeWitt 


foils'/ 


r  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Surface  Weapons  Center  (Code  K51) 
Dahlgren,  Virginia  22448 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  ft  WORK  UNIT  NUMBERS 


OMN : 0K50TT00 1 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS  f  RE  PORT'D  ATE - - -  . 

Strategic  Systems  Projects  Office  (  II  }  JanR)MlEB8y  ( 

Washington,  DC  20390  . 

14.  MONITORING  AGENCY  NAME  A  ADDRESS*1//  dlllarmnl  Tram  Controlling  Offlepj  IS.  SECURITY  CLASS,  (a!  I Ma  raport) 


UNCLASSIFIED 

IS*.  DECL  ASSI  FI  CATIOn/ DOWNGRADING 
SCHEDULE 


I  IS.  DISTRIBUTION  STATEMENT  (at  Ihl.  Raport) 


Approved  for  public  release;  distribution  unlimited. 


I  17.  DISTRIBUTION  STATEMENT  (ol  lha  abstract  antarad  In  Block  70,  II  dllloront  from  Raport) 


I  IS  SUPPLEMENTARY  NOTES 


JUf  r  y 

Cr<  \'fo  "  r  ' 

V.  v 

V  a  •  r 


1 19.  KEY  WORDS  (Continue  on  rmrmram  mid*  1/ nictMify  and  Identity  by  block  ftuotbtO 


20.  ABST^CT  (Contlnum  on  rmrarma  ml  dm  It  nacammmry  *nd  Identity  by  block  mambmt) 

The  equation  of  motion  for  the  restricted  two-body  problem  is  trans¬ 
formed  via  the  Kustaanheimo-Stiefel  transformation  method  (KS)  into  a 
dynamical  equation  in  the  standard  four-space  representation.  Pure 
rotations  about  an  arbitrary  axis  are  introduced  using  the  Euler  parameters 
which  are  then  generalized  to  include  dilations.  Kinematic  equations  in 
the  generalized  Euler  representation  are  compared  to  the  dynamical  equations 
of  the  KS  representation  giving  a  simple  interpretation  of  the  latter* 


FORM 
l  JAN  73 


EDITION  OF  1  NOV  SS  IS  OBSOLETE 

S/N  0!  02- LF -01 4-6601 


UNCLASSIFIED  j 

SECURITY  CLASSIFICATION  OF  THIS  PAO« 


SECURITY  CLASSIFICATION  OF  this  PAGE  (Whit  Dim  Bnfr.d) 


FOREWORD 


The  work  represented  by  this  technical  report  was  performed  in  the  Fire 
Control  Presetting  Analysis  Branch  of  the  FBM  Geoballistics  Division,  Strategic 
Systems  Department  under  Task  Assignment  No.  37430. 

This  report  has  been  reviewed  by  J.  R.  Fallin  and  D.  L.  Owen  of  the  FBM 
Geoballistics  Division. 


R.  T.  RYMND  JR.',  Head 
Strategic  Systems  Department 


111 


hintfffii  *»  ''**  ;i *-*•*■• 


ACKNOWLEDGEMENT 


Mr. 


The  author  wishes  to  express  his  appreciation  to  Dr.  J.  Blanton  and 
0.  Schultz  £or  reading  parts  of  the  manuscript. 


IV 


INTRODUCTION  . 

EQUATION  OF  MOTION  . 

KS  TRANSFORMATION  MATRIX  .... 
ROTATION  ABOUT  AN  ARBITRARY  AXIS 

CONCLUSION  . 

REFERENCES  . 

DISTRIBUTION 


INTRODUCTION 


Stiefel  and  Scheifele1  have  given  a  general  discussion  of  the  Kustaanheimo- 
Stiefel2  transformation  method  (KS)  in  the  two-body  problem.  Many  papers  have 
appeared  in  which  specific  problems  or  applications  have  been  addressed;  but  the 
feeling  persists  that  in  going  from  a  three-space  to  a  four-space  representation, 
fundamental  understanding  is  impaired.  Thus,  there  seems  to  be  a  reluctance  to 
apply  this  method  to  trajectory  determination,  guidance  algorithms,  or  rigid-body 
kinematics  in  spite  of  favorable  claims  by  some  authors  with  regard  to  its  effi¬ 
ciency  and  accuracy.3*4'5 

It  is  the  objective  of  this  paper  to  show  that  the  KS  method  can  be  under¬ 
stood  in  terms  of  a  generalized  Euler  parameterization  of  rotations  and  dilations. 
The  second  section,  EQUATION  OF  MOTION,  sets  up  the  equation  of  motion  for  the 
restricted  two-body  problem  while  the  third  section,  KS  TRANSFORMATION  MATRIX, 
briefly  reviews  the  KS  method  and  applies  it  to  the  two-body  equation  of  motion. 
The  fourth  section,  ROTATION  ABOUT  AN  ARBITRARY  AXIS,  sets  up  the  Euler  parameter 
description  of  rotations  and  then  generalizes  it  to  include  dilations.  Finally, 
the  kinematical  equations  of  this  representation  are  compared  to  the  KS  dynamical 
equations  of  motion  to  show  that  they  have  equivalent  interpretations. 


EQUATION  OF  MOTION 


The  equation  of  motion  per  unit  mass  in  the  restricted  case  is 

r  +  M  r  =  F  (r,  r,  t)  (1) 

where  p  =  GM  and  F  is  the  perturbing  acceleration.  If  E  is  the  eccentric  anomaly, 
then  consider  a  scaled  eccentric  anomaly  defined  by 


s 


(2) 


1 


where  V  is  the  geometric  mean  of  the  Keplerian  velocities  at  apocenter  and  peri- 

ID 

center.  Then  a  change  of  variables  can  be  made  with 

dt  =  rds  (3) 

so  that  the  equations  of  motion  become 

“2  ° r5  ”  ~3  ?  +  r  =  F  (r,  r,  s)  (A) 

where  r  denotes  the  derivative  with  respect  to  s.  In  the  case  where  F=0  (the 

<v 

Kepler  problem),  the  energy  is  a  constant  given  by 

Ek  =  Ekinetic  +  potential  =  "  r  =  "20fk 

where  >  0  for  elliptical  orbits,  the  only  case  of  interest  here. 


KS  TRANSFORMATION  MATRIX 

P.  Kustaanheimo  and  E.  Stiefel2  proposed  a  regularization  method  by  intro¬ 
ducing  a  A  x  A  transformation  matrix  and  four-component  vector  given  by 


L(u)  = 


U1 

•U2 

'U3 

u, 

U2 

U1 

'UA 

-a, 

u3 

UA 

U1 

u, 

UA 

"U3 

U2 

-u. 

such  that  r  =  L(u)u;  that  is 


/ 


r  = 


3 
0 

-  ,.2 


=  L(u)u  = 


2  2  2  2 

VWU4 

2(»jU2  -  u3u4) 
+  u2u4) 


\ 


/ 


(7) 


and  |r  |  =  u  *  u  =  uz.  Since  L(u)  is  an  orthogonal  matrix,  the  inverse  L  *(u)  is 


related  to  the  transpose  Li(u)  by 


L-1(u)  =  LT(u) 

^  «x» 


(8) 


Thus 


o  1  tT,  \  o 
u  ~  2u*  L  (u)  r 

<V  IV  M 


(9) 


with  u  and  u  satisfying  a  bilinear  relation,  denoted  £(u,  u)  =  0,  such  that 


ACu,  u)  =  u^  -  u3a2  +  u2S3  -  UjU4  =  0 


(10) 


and  £(u,  u)  is  a  constant  of  the  motion  in  the  perturbed  two-body  problem.  The 

details  of  this  method  have  been  discussed  by  E.  Stiefel  and  6.  Scheifele.1 

Introduction  of  the  KS  transformation,  equation  (7),  into  the  equations  of 
motion  (4)  leads  to  the  u-space  equation  of  motion 


”*  <  h  ■  r  > ;  ■ §  iT(">  I 


(in 


which  with  the  energy  relation  (S)  becomes 


oo 


u  =  |  LT(u)  F 


(12) 


Recall  that  in  the  Keplerian  problem  the  potential  energy  is 


E  =  V  =  -  y 

potential  kep  r 


so  that  the  total  Keplerian  energy  is 


E  =E  +  v  =  —  -  ^=  -2o( 
Kep  Kin.  Vkep  2  r  ^"k 


The  total  energy  Efc  can  be  written  as 

v2  v2 

E  =  E,  .  +  E  =  x-  +  V,  +  V  „  +  V  =  ^-  +  V 

t  kin.  pot.  2  kep  pG  p  2  t 


where 


V  =  '  v/t 

VpG  =  Potential  due  to  perturbing  (non-spherical)  gravity 

Vp  =  Potential  due  to  non- gravitational  perturbations 

V.  =  Total  Potential  =  V,  +  V  „  +  V 
t  kep  pG  p 


Then  Efc  can  be  written  as 


=  -2a.  +  V  _  +  V 
t  k  pG  p 


With  the  definition  2a^  =  -E^  the  above  becomes 


a  ~  a  -  -  V  -  -  V 
0  k  2  pG  2  p 


(13) 
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Now  the  total  perturbing  acceleration  F  can  be  split  into  that  derivable  from  a 
potential  V  plus  that  which  is  not  (P) ,  so 

F  =  -  3  V  +  P 
~  r  tP  ~ 

—  M  - 

=  -a  v  .  -  a  v  +  p 

r  pG  r  P  ~ 


The  equations  of  motion  (12)  then  become 


or 


“  *  V  = '  I V I '  I  vp  t*  i  iT(u)  i 


°u  +  «-u  =  -  |  V  _  u  -  5  V  u  +  Q 


where 


Q  =  f  lt(u)  f  = 

>w  ^  #v  «w 


I  8rvtP  + !  iT<»>  i 


It  can  be  seen  that 


LT(u)  a  =  |  3  V 
-  r  2  u 


by  expanding  both  sides,  so  that  the  equations  of  motion  become 


?  ♦  “or  •  I  v r  I Y *  § (-  K  v  ♦  r 


(14) 
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If  all  perturbations  are  derivable  from  a  potential  V  ,  then  P  =  0  and  no  terms 
with  L  (u)  appear  in  the  equations  of  motion  giving 


°u  +  a.u  = 

<v  ''/v 


-Iv.u-rV  u  -  7  a  v  « • f  a  v 
2  pG  ~  2  p  ^  4  u  pG  4up 


(15) 


If  all  perturbations  are  due  to  the  gravitational  field,  then  =  0  and 


u+  o.u  =  ■  J  V  ,  u  •  7  a  V  _ 
~  2  pfl  „  4  u  pG 


(16) 


du 


Finally,  some  authors  use  u*  =  -r=r  so  that  with  ur  =  a 

at  o  o 


u  =  2u>  u* 
o 


and  equation  (16)  becomes 


u**  +  i  u  =  77—9 — (V  „  u  +  7  3  V  „) 
a  4  “  pG  j,  2  u  pG'' 


(17) 


Using  equation  (12)  for  the  equation  of  motion,  a  substitution  Of^  =  u>£  and 


r  T 

j  L  (u)F  =  D(u)  gives  the  form 


°i?  +  u>£  u  =  D(u) 

**  ^  *v  ^ 


which  is  the  form  for  driven,  coupled  oscillators. 
D(u)  =  0,  the  homogeneous  solution  is 


(18) 


'0) 

If  u'  is  a  solution  for 


u^(s)  =  u^(0)  cos  (w.s)  +  —  u^  (0)  sin  (u».s) 


(19) 
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while  the  general  solution  is  the  integral  equation 


u(s)  =  u^(s)  +  -j-  sin  (u»  s) 


(u».o)  D  (u(o))do 


(20) 


1_ 

U) 


cos 


<v> 


i 


sin  (u)ka) 


D  (u(o))do 


A  method  of  solution  is  to  use  the  unperturbed  solution  u^^(s)  in  the  integrands, 

the  first  Born  approximation,  leading  to  a  set  of  elementary  integrals  to  evaluate. 
Andrus3  has  done  this  with  D  for  the  gravitational  perturbation  acceleration 
including  oblateness. 


ROTATION  ABOUT  AN  ARBITRARY  AXIS 

If  a  vector  r  is  rotated  positively  about  an  axis  n  thru  an  angle  0,  then  the 
resulting  vector  r'  is  given  by 

r'  =  r  cos  0  +  (1  -  cos  8)(n  •  r)  n  -  sin  0  (r  x  n)  (21) 

This  is  equivalent  to  a  negative  rotation  of  the  axis.  The  former  rotation  is 
referred  to  as  an  active  rotation  while  the  latter  a  passive  rotation.  If  a 
coordinate  system  is  introduced  such  that 

n  =  (cos  a,  cos  p,  cos  y) 

r  —  (Xj,  ,  x^) 


7 


then  equation  (21)  becomes 


r'  =  (I  -  2  sin2  |  A  (a,p,Y)  +  2  sin  |  cos  f  B  (ot,P,Y)]r 


where 


A(a,B,Y)  = 


sin2a 


cos  a  cos  $ 


cos  a  cos  Y 


-  cos  oi  cos  P 
sin2p 

-cos  P  cos  Y 


-  cos  a  cos  Y 


-  cos  P  cos  Y 
8in2Y 


B(a,p,Y)  = 


cos  Y 
-  COS  P 


-  cos  Y 


cos  a 


cos  p 


-  cos  a 


and  I  is  the  unit  matrix 


Equation  (22)  is  the  Whittaker6  result.  Introducing 


the  Euler  parameters 


qo  =  cos  j 

.  e 

q^  =  cos  a  sin  ^ 

„  •  6 
q^  =  cos  p  sin  n 

6 

q3  =  cos  Y  810  2 


with  2  q?  =  1,  the  above  becomes 
A 


qf  "  q2  ”  q3  +  q0 


r’  =  I  2(qjq2  +  q3qQ) 

\2(q1q3  - 


2(qjq2  *  Q3V 
-q2  +  q|  "  ^3  +  q0 
2(q2q3  +  q^Q) 


2(q1q3  +  q2V  \ 

2(q2q3  -  qjqo)  J  l 
.  q2  -  q2  +  q|  ♦  qg  /  (25a) 


f 


or 


r'  =  R(q)  r  (25b) 

^  M 

sw 

A  negative  rotation  of  the  vector  (or  positive  rotation  of  axes)  produces  the 
transpose  of  the  matrix.  This  can  be  seen  directly  since  pairs  of  q's  in  the  off 
diagonal  terms  which  do  not  contain  qQ  are  even  in  0  while  those  containing  qQ 
are  odd. 

A  triad  of  orthogonal  vectors  can  be  introduced  by 


x  = 


fWq3+q0 
2(qiq2-q3q0)  |  .  y  = 

,2(qlq3+q2q0^ 


/  2(qlq2+q3q0}> 
~ql+q2’q3+q0 

2(q2q3+qiq0) 


z  = 


^2(qlq3'q2q0^ 

2(q2q3+q1q0) 

-qf-q|+q!+qgi 


(26) 


Thus  equation  (25)  represents  a  projection  of  r  onto  the  rotated  axes  x,  y,  z  and 


the  rotation  is  generated  by 

*  =*£(q)q>  y  =/(q)p.  z  =jt(q)s 

<v  m  m  ^  ^  ^ 


(27) 


(28) 


(29) 
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Note  the  identity  between  the  x  vector  and  the  position  vector  r  =  L(u)u, 

£*  M  M  M 

Me 

equation  (7),  of  the  KS  transformation  if  and  q^  =  u^,  i  =  1,2,3.  One 

could  introduce  two  other  vectors  orthogonal  to  r  in  the  KS  formalism  as 


2(ulVu3U4) 

U4~U1+U2~U3 

2(U2U3’U1U4) 


2(U1U3_U2U4) 

2^U2U3+Uiu4) 

u2-u;-u2+u| 


and  generated  by 


r_.  =  L(u) 

fS0  ^  /v 


,  rM  =  L(u) 


and  interpret  the  KS  formalism  as  a  parameterization  of  a  rotated  set  of  axes  onto 
which  an  initial  vector  is  projected  to  obtain  the  resulting  vector. 

For  completeness,  a  fourth  vector 


2d  = 


can  be  introduced  which  is  orthogonal  to  q  ,  q.  and  q  .  These  orthogonal  vectori 

»a  .vC 

4  =  -  =  .  , 

now  span  R  .  In  particular,  it  is  not  hard  to  show  that  q  can  be  expanded  with 

a*  / 


% 


components  along  q.  ,  q  and  q.  which  results  in 

~d 

i  =  i  !M  3. 


where 


B(u»)  = 


and  ut  =  (u»j,  u^,  u^)  is  the  angular  velocity  vector  and  uk  =  0^.  Note  that 
equation  (33)  can  also  be  written  as 


|a  =  2 


From  equation  (33)  it  is  seen  that  with  q  <->  q 
q  =  5  (B  q  +  B  q)  =  !  (B  q  +  i  B2  q) 


But  with  the  fact  that  B2  —  -  u>2  I ,  there  results 


ujz  i  • 
q+/— q  =  ^Bq 


f 


Recall  the  change  of  variable  introduced  in  equation  (3)  giving  for  a  particle  a 
distance  r  from  the  origin 


If,  for  the  time  being,  r  is  considered  constant  for  pure  rotations  of  the  coor¬ 
dinate  system  or  circular  motion  of  the  particle,  equation  (36)  becomes 


oo 

q  + 


r2w2 


A 


(37) 


But  for  pure  rotations,  the  total  energy  is  composed  of  the  rotational  energy 
^  r2u)2  and  the  potential  energy  -p/r,  so  there  results 

Ekep  =  ‘2ak  =  \  r2“2  -  M/f  *  5  r2«>2  -  r2  u»2  =  -  u>2 


or 


r2u»2 

A 


and  equation  (37)  becomes 


oo 

q  + 


R  0* 


oo_  r  ' 
q  -  -  b 


2  "  q 
<v  IV 


(38) 


which  has  the  same  form  as  equation  (12).  The  case  where  r  f  0  will  be  addressed 

later.  However,  at  this  point,  it  is  suspected  that  the  KS  method  is  similar  to 

o 

the  Euler  parameterization  of  rotations,  and  that  B  is  related  to  the  perturba- 
tions.  Z 
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Another  representation  of  equation  (25)  is  obtained  by  introducing  the  qua¬ 
ternion  notation  such  that 

0  s  %  *  9lel  *  92e2  *  ’3e3  =  «o  +  | 

r  =  0  +  Xjej  +  *2e2  +  x3e3  =  0  +  *  (39) 

Q"1  =  ^  ’  l 

Using  quaternion  multiplication,  equation  (25)  becomes 

r'=Qr  Q1  (40) 


Note  that  n  is  an  eigenvector  of  the  rotation  matrix  R  (q) ;  that  is 
R(q)  ft  =  ft 


(41) 


or  in  quaternion  representation 
QftQ  1  =  ft 


(42) 


Recall  that 

I  qj  =  1  X  =  0,  1,  2,  3 

so  that  the  transformation  represents  pure  rotations.  This  can  be  generalized7 
by  introducing 

u\ =  5  q\  (43) 


13 


/ 


with  u2=2u£  =  62.  Then  introducing  the  quaternion 

U  =  uQ  +  u  =  6^  +  q)  =  6Q  (44) 

U_1  =  0*  (u0  *  f 

there  results 

r'  =  UrU-1  (45) 

and  the  transformations  now  contain  rotations  and  dilations.  Using  the  matrix 
representation,  equation  (22)  can  be  written  as 

r'  =  R(u)  r  (46) 

a*  m*  A# 


where  R  (u)  is 

A# 


the  matrix  in  equation  (25)  with  q^  =  u^  and  without  the  restric¬ 


tion  that  2  uj*  =  1.  Also,  if  u^  is  associated  with  qQ  in  the  x  vector  in  equa¬ 
tion  (26),  there  results  the  r  vector,  equation  (7),  of  the  KS  method.  Thus  it 

may  be  considered  that  the  KS  method  is,  in  essence,  an  Euler  parameterization  of 
transformation,  generalized  in  the  sense  of  Velte.7  The  transformations  are  gen¬ 
erated  in  the  KS  method  by 

r  =  l(u)  u  , 

A#  #V  M 


in  this  paper  by 

x  =  «£(u)  u 

A/  A/  A# 
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(47) 


or  as  in  Vitins8  by 
x  =  A(u)  u 


where 


A(u)  = 


u  = 


h 

'  U2 

V4 

“A 

U2 

"3 

V 


”U2 

'U3 

u 

i 

U1 

u4 

u. 

3 

1 

U1 

-u, 

U3 

'U2 

-u 

(48) 


The  equation  of  notion  for  q  (pure  rotation)  is  given  by  equation  (37).  With 

g 

u.  =  6q.  ,  the  equation  of  motion  for  u  (which  now  includes  rotations  and  dila- 
/v  A  m 

tions)  becomes 


oo  y  2|||2  OO  O  t»  O 

u  +  — t —  u  =  [6/6  I  +  2  66  IB  (u>)  +  x  B  (u>)]  u 

<w  18  <w  ***  w  •  »v  M 


(49) 


For  pure  rotations  6  =  constant  and  equation  (49)  reduces  to 


oo  .  r2u>2  _  r  8,  ,, 

u  +  — t  u  =  x  B(u>)  u 

«W  **  “  »W  M 


or 


oo  r  * 

q  +  — T~  q  =  f  b(w)  q 
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which  is  equation  (37).  Using  the  relation 


where 


r2ui2 

~7T  =  °k  ’  6 


e  =  r  (J*  - 

c  =  area  constant 
b  =  semiminor  axis 
equation  (49)  becomes 


nn  OO  °  v  O 

u  +  a^u  =  [6/6  1  +  2  66B(uj)  +  ^  B(w)  +  A)  u 


(50) 


Comparing  equation  (50)  with  equation  (12)  it  is  seen  that  the  kinematical  equa¬ 
tion  due  to  the  Euler  parameterization  of  rotations  and  dilations  has  the  same 
form  as  the  dynamical  equation  of  the  KS  method. 

CONCLUSION 

A  simple  interpretation  of  the  KS  method  can  be  made;  namely,  the  method  is 
equivalent  to  treating  the  dynamics  in  a  rotating  and  dilating  coordinate  system 
characterized  by  generalized  Euler  parameters. 
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